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ABSTRACT 

We present high-sensitivity 2'x4' maps of the J=2— >1 rotational lines of SiO, 
CO, 13 CO and C 18 0, observed toward the filamentary Infrared Dark Cloud (IRDC) 
G035. 39-00. 33. Single-pointing spectra of the SiO J=2— >l and J=3— >2 lines toward 
several regions in the filament, are also reported. The SiO images reveal that SiO 
is widespread along the IRDC (size ^2pc), showing two different components: one 
bright and compact arising from three condensations (N, E and S), and the other 
weak and extended along the filament. While the first component shows broad lines 
(linewidths of ^4-7 km s -1 ) in both SiO J=2— >1 and SiO J=3— >2, the second one is 
only detected in SiO J=2— >1 and has narrow lines (~0.8 km s -1 ). The maps of CO 
and its isotopologues show that low-density filaments are intersecting the IRDC and 
appear to merge toward the densest portion of the cloud. This resembles the molecular 
structures predicted by flow-driven, shock-induced and magnetically-regulated cloud 
formation models. As in outflows associated with low-mass star formation, the exci- 
tation temperatures and fractional abundances of SiO toward N, E and S, increase 
with velocity from ^6 to 40 K, and from ^10~ 10 to ^10~ 8 respectively, over a velocity 
range of ~7kms _1 . Since 8^m sources, 24 /zm sources and/or extended 4.5 /mi emis- 
sion are detected in N, E and S, broad SiO is likely produced in outflows associated 
with high-mass protostars. The excitation temperatures and fractional abundances of 
the narrow SiO lines, however, are very low (^9K and ^10~ u , respectively), and 
consistent with the processing of interstellar grains by the passage of a shock with 
u s ^12kms _1 . This emission could be generated i) by a large-scale shock, perhaps 
remnant of the IRDC formation process; ii) by decelerated or recently processed gas 
in large-scale outflows driven by 8 fim and 24 jirsx sources; or iii) by an undetected and 
widespread population of lower mass protostars. High-angular resolution observations 
are needed to disentangle between these three scenarios. 
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1 INTRODUCTION 

Infrared Dark Clouds (IRDCs) are high-extinction regions 
viewed against the diffuse mid-IR Galactic background 
JPerault et al.lll996l;lEgan et aljl998h . These clouds are cold 
(T <25K: |Pillai et al.ll2007l ) and exhibit a range of densities 
from n (H)^10 3 cm" 3 to ^10 4 -10 5 cm 3 in their clumps and 
cores jTevssier. Hennebelle fc Peraultl 120021 : iButler fc Tan] 
2009). Since these structures have masses and mass surface 
densities similar to regions that are known to be forming 
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massive protostars and star clusters, they may represent the 
initial conditions for massive star and star cluster forma- 
tion (Rathborne, Jackson fc Simon 2006; IZhang et al.ll2009l ; 
iRagan. Bergin fc Gutermutbll2009l h 

It is well-known that silicon monoxide (SiO) is an ex- 
cellent tracer of molecular gas processed by shocks. While 
SiO is heavily depleted onto dust grains in the quiescent gas 
of dark clouds such as L183 (upper limits of the SiO frac- 
tional abundance of ^10~ 12 ; Ziurys, Friberg & Irvine 1989; 
iRequena- Torres et al.l [20071 ). this molecule is enhanced by 
large factors (in some cases by >1 6 ) toward molecular out- 
flows (Ma rtin-Pintado et al.|[l992r i. This is due to the injec- 
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Table 1. Observed molecular transitions and line frequencies, 
telescope beam sizes and beam efficiencies for the EMIR and 
HERA receivers at the IRAM 30 m telescope. 



Transition 


Frequency (MHz) 


Beam size (") 


Beam Eff. 


SiO J=2^1 


86846.96 


28 


0.81 


SiO J=3^2 


130268.61 


19 


0.74 


CO J=2->1 


230538.00 


11 


0.63 


13 CO J=2->1 


220398.68 


11 


0.52 


C ls O J=2->1 


219560.36 


11 


0.52 



tion of molecular materi al into the gas phase by the process- 
ing of dust grains (e.g. Caselli. Hartquist fc Havnesl 19971; 
Schilke et all ll997: Jimenez-Serra et al.ll2008l : iGuillet et al l 



20071. 120091) . 

The typical SiO emission measured in molecular out- 
flows show s broad line profiles with lin ewidths of some tens 
of kms -1 (|Martm-Pintado et ai1ll992T ). Narrower SiO lines 
have also been detected toward low-m ass star forming re- 
gions such as NGC1333 a nd L1448-mm (|Lefloch et al.ll 19981 ; 
I Jimenez-Serra et all2004 | ). Although the nature of this emis- 
sion is not clear vet. lLefloch et al"] j 19981 ) have proposed that 
these lines could trace shocked material deflected and decel- 
erated by the interaction with pre-existing clumps. Alterna- 
tively, narrow SiO could arise from gas recently processed 
by the magneti c precursor of young magne to- hydrodynamic 
(MHD) shocks (| Jimenez-Serra et al-lliool). 

In the case of IRD Cs, iMotte et ail (|2007l ) and 
iBeuther fc Sridharanl (|2007h have recently carried out two 
large surveys of SiO emission toward the Cygnus X molec- 
ular cloud complex and toward a sample of IRDCs, respec- 
tively. Their single-pointing observations show that the de- 
tection rate of SiO toward IR-quiet massive cores close to 
sites of on-going star formation, is relatively high. This is 
expected since SiO is tightly associated and restricted to 
shocked gas in outflows. However, no widespread narrow 
SiO emission, as seen toward NGC1333, has been reported 
in IRDCs so far. 

We present the detection of widespread SiO J=2— s>l 
emission (si ze of ^2pc) towar d the IRDC G035. 39-00. 33 
(Cloud H in lButler fc Tanlfeoogl ). From its filamentary mor- 
phology, this IRDC is believed to be at the early stages of its 
evolution as predicte d by dynamical models of giant molec- 
ular cloud formation (|van Loo et al.ll2007l ; iHcnncbcll e et al.l 
|200S| , Heitsch, Stone & Hartmann 2009). The observed large- 
scale SiO feature is probably a composition of broad and 
compact emission, linked to outflows associated with high- 
mass star formation; and extended narrow SiO emission. 
The observed narrow SiO lines could be explained i) by a 
large-scale shock, remnant of the IRDC formation process; 
ii) by decelerated or recently processed gas in the precursor 
of MHD shocks in large-scale outflows, probably driven by 
the 8 /im and 24 /im sources observed in the IRDC; or iii) 
by an undetected and widespread population of lower-mass 
protostars. 



2 OBSERVATIONS AND RESULTS 

The J=2^1 lines of SiO, CO, 13 CO and C 18 0, were 
mapped with the IRAM (Instituto de Radioastronorma 
Milimetrica) 30 m telescope at Pico Veleta (Spain) over 
an area of 2'x4' toward G035. 39-00.33. These observa- 
tions were carried out in August 2008, and in January and 
February 2009. The large-scale molecular images were ob- 
tained in the On-The-Fly (OTF) mode using the offsets 
(1830",658") for SiO, 13 CO and C 18 0, and (4995",2828") 
for CO, as off-positions. The central coordinates of 
the map were a(J2000)=18 ,I 57 m 08 s , <5(J2000)=02° 10'30" 
(1=35.517°, b= -0.274°). The SiO J=2->1 emission was 
mapped with the old SIS receivers, while the HERA multi- 
beam receivers simultaneously observed the J=2— >1 tran- 
sitions of 13 CO and C ls O. The CO J=2^1 emission was 
mapped with the new generation EMIR receivers. In addi- 
tion, we carried out single-pointing observations of the SiO 
J=2— >1 and J— 3— 5-2 emission with EMIR toward the off- 
sets (0,20), (30,-30), (0,-80), (-10,70) and (-37,37). The for- 
mer three positions correspond to the brightest SiO emission 
peaks observed toward the IRDC (Section [2~Tj) . The latter 
two offsets show the regions where we have detected narrow 
SiO lines (Section l2.2p . All receivers were tuned to single 
sideband (SSB) with rejections of ^10 dB. The beam sizes 
were 28" at 90 GHz for the SiO J=2^1 line, 19" at 130 GHz 
for the SiO J=3^2 emission, and 11" at 230 GHz for the 
CO, 13 CO and C 18 J=2^1 lines. The VESPA spectrom- 
eter provided spectral resolutions of 40 and 80 kHz, which 
correspond to velocity resolutions of ~0.14 and 0.1 kms - 
at 90 and 230 GHz, respectively. Saturn was used to calcu- 
late the focus, and pointing was checked every two hours on 
G34.3+0.2. Typical system temperatures ranged from 100 
to 300 K. All intensities were calibrated in units of antenna 
temperature, T%- To convert these intensities into units of 
main-beam temperature, T m t, we have used efficiencies of 
0.81, 0.74 and 0.63 for the EMIR data at ~90, 130 and 
230 GHz, and of 0.52 for the HERA data at ~230 GHz. All 
this information is summarized in Table[T] 

Figure]]] (central panel) presents the high-sensitivity 
map of the SiO J=2— >1 emission integrated from 36 to 
51.5 kms -1 (blue contours), and superimposed on the mass 
surface dens ity map, with an a ngular resolution of 2", 
reported by iButler fc Tanl (|2009h ■ The SiO J=2-»l map 
has been obtained by averaging the OTF dumps in the 
SiO beam of 28", and by using a Nyquist-sampled grid 
with a pixel size of 14". The 2cr intensity level of SiO is 
shown in dashed contours. The l ocation of the massi ve cores 
(crosses; [Rathborne et al]|2006] ; [Butler fc Tanll2009l ). 24/jm 
sour ces (red open trian gles; extracted from MIPSGAL im- 
ages; ICarev et alj|2009h . 8 /im sources (red open circles), and 
4. 5 /im extended emission (green squares; called green fuzzies 
in lChambers et al. I l2009l ) in this IRDC, are also shown. The 
flux lower limits of the 8 ^m and 24 /im sources reported 
in Figure]]] are >3.5mJy and ^2mJy, respectively. We note 
that the cavity- like structures seen around the 8 /im sources, 
are pro duced by the fact tha t the extinction mapping tech- 
nique of lButler fc Tanl l|2009l) cannot be applied in the vicin- 
ity of IR-bright sources. 

From Figure[T] we find that the SiO J=2—¥l emission 
is widespread across the filament with a spatial extent of 
>150"x50". This corresponds to >2.1pcx0.7pc at a dis- 
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Figure 1. Central panel: Integrated intensity map of the SiO J=2— ^1 lin e toward G035 . 39-00. 33 for the velocity range from 36 to 
51.5kms — 1 (blue contours), overlapped on the mass surface density map of lButler fe Tar] | |2009| , gray scale). The contour levels of the 
SiO emission, are 10 (2a; dashed contour), 15, 20, 30 and fflrnKkms" 1 . For the mass surface density map, contours are 0.014 (2cr) , 
0.02 1, 0.035, 0.049 , 0.07, 0.105 and 0.14gcm — 2 , respectively. Crosses indicate the cores reported in the cloud bv lRathborne et al . (2006) 
and iButler fc Tanj J2009I) . Red open circles, red o pen triangles and green squares show the location of the 8/im sources, 24 /an sources 
l| Carev eT*aTTl2*009l ) and 4.5 /an extended emission ijChambers et ai]|200Sjh detected in G035. 39-00. 33, respectively. The marker sizes used 
for the 8 /im and 24 /im sources, have been scaled by the source flux. Red dashed lines show the directions of the P — V diagrams of 
Figure[2] Red numbers indicate the distance (in arcseconds) along the cuts made for these diagrams. The SiO J=2— >l beam size is plotted 
at the upper right corner. Left and right panels: Sample of spectra of SiO J=2— >1 and J=3— >2, and of CO, 13 CO and C 18 J=2— >1, 
measured toward several positions across the IRDC. The CO, 13 CO, and C 18 spectra were obtained by averaging the OTF data within 
the 28"-beam of the SiO J=2— >1 observations. The vertical dotted lines show the velocities of the CO filaments, as calculated from the 
C ls O J=2-¥l emission (Tableland SectionEOl) . 



tance of ~2.9kpc ()Rathborne et al.ll2006h . This size should 
be considered as a lower limit since SiO is also detected 
toward the north and northwest of Core MM7/H6 (see 
Section [2~2"1l , and extends off to the south of the imaged area. 
The SiO map shows three bright condensations (marked N, 
E, S) surrounded by weaker, more extended emission, that 
covers much of the IRDC filament. 

In FigureQ] we also report a sample of spectra of SiO 
J=2^1 and J=3->2, and of the CO, 13 CO and C 18 
J— 2— >1 lines, measured toward several positions in the 
IRDC. Since single-pointing observations of SiO J=2— 5-1 
and J=3— >2 were not carried out toward (10,10), the SiO 
J=2—Yl spectrum reported in Figure[T]was obtained by av- 
eraging the OTF data within the 28"-beam of the IRAM 
30 m telescope at ~90GHz. The spectra of CO and of its 
isotopologues shown in Figure[T] have also been obtained fol- 
lowing this procedure. From the SiO spectra, we find that 
SiO shows a wide variety of line profiles, from broad emission 



with red- and/or blue-shifted line wings [see offsets (0,20), 
(30,-30) and (0,-80)], to narrow lines peaking at the ambient 
cloud velocity tJ iSfl ~45kms _1 [offsets (10,10), (-10,70) and 
(-37,37)]. The CO, 13 CO and C ls O lines show three different 
velocity components centred at WLSfl~44, 45 and 47 kms -1 
(see Section r2~31 . 

In Table[2] we report the observed parameters (cen- 
tral radial velocity, v^sr, linewidth, Av, and peak inten- 
sity, TjJ) of the different velocity components measured in 
SiO J=2->-l, SiO J=3^2 and C ls O J=2^1 (representa- 
tive of the low-density CO gas) toward (0,20), (10,10), (30,- 
30), (-10,70), (-37,37) and (0,-80). These parameters were 
obtained by fitting the molecular line emission with three 
gaussian line profiles simultaneously. Although this method 
works better for the C ls O J—2~ >1 emission, the gaussian 
linewidths derived from the SiO 3=2— >\ line profiles pro- 
vide a rough estimate of the velocity range extent and ter- 
minal velocities of the shocked SiO emission. In Table(2] the 
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Table 2. Observed parameters of the SiO 3=2— >l and J=3— >2 lines, and of the 3=2— >1 emission of C 18 0, 
toward several offsets in IRDC G039.35-00.33. 



Molecule OFFSETS 







(0,20) 






(10,10) 






(30,-30) 






VLSR 


Av 


T*„ 
,1 




Av 


T*„ 

A 


U LSR 


Av 


T*. 
A 




(km s — ) 


(kms -1 ) 


(K) 


(kms -1 ) 


(km s — ) 


(K) 

V- 1 - V 


(kms -1 ) 


(kms -1 ) 


(K) 


SiO(2->l) 


-44.0 




sJO.03 


44.73(8) 


0.7(2) 


0.07(1) 


42.9(4) 


7.0(8) 


0.06(1) 




45.56(6) 


1.3(2) 


0.13(1) 


45.7(1) 


0.8(3) 


0.06(1) 


44.73(7) 


1.9(2) 


0.12(1) 




47.30(8) 


1.0(2) 


0.07(1) 


-46.8(2) 




<0.04 


-46.4 




sSO.04 


SiO(3^2) 


-43 




<0.04 








42.9(3) 


3(1) 


0.06(2) 




45.5(1) 


1.7(4) 


0.07(1) 








44.9(2) 


2.3(6) 


0.08(2) 




47.0(2) 


0.8(3) 


0.05(1) 








—46.4 




^0.05 


C 18 0(2^1) 


44.0(2) 


2.3(3) 


0.34(3) 


44.500(1) 


3.0(4) 


0.25(4) 


43.1(2) 


1.7(4) 


0.08(3) 




45.41(1) 


1.08(4) 


1.30(3) 


45.38(1) 


1.07(3) 


1.36(4) 


45.21(2) 


1.30(3) 


0.87(3) 




46.77(2) 


1.26(5) 


0.78(3) 


46.78(8) 


1.4(1) 


0.36(4) 


46.4(1) 


1.4(3) 


0.09(3) 






(-10,70) 






(-37,37) 






(0,-80) 




SiO(2-+l) 


-44.9 




sSO.018 


-43.9 




^0.03 


43.2(3) 


6.2(3) 


0.04(1) 




-45.9 




^0.018 


45.3(1) 


0.8(3) 


0.04(1) 


45.6(3) 


2.9(3) 


0.11(1) 




46.5(2) 


1.5(4) 


0.021(6) 


-46.7 




^0.03 


49.1(3) 


3.5(3) 


0.06(1) 


SiO(3->2) 


-44.9 




<0.04 


-43.9 




^0.05 


-43.2 




sS0.07 




-45.9 




<0.04 


~45.1 




^0.05 


45.3(1) 


0.9(3) 


0.07(2) 




-46.9 




^0.04 


-46.7 




^0.05 


47.5(5) 


4.8(9) 


0.11(2) 


C 18 0(2^1) 


44.9(1) 


2.2(2) 


0.57(3) 


43.9(1) 


1.1(1) 


0.69(3) 


44.097(8) 


2.4(4) 


0.13(3) 




45.89(3) 


0.8(1) 


0.60(3) 


45.1(1) 


1.2(1) 


0.33(3) 


45.07(2) 


1.22(3) 


1.07(3) 




46.92(5) 


1.34(9) 


0.68(3) 


46.7(1) 


1.7(1) 


0.40(3) 


45.92 (3) 


1.04(6) 


0.56(3) 



errors in vlsr and Av correspond to those obtained from 
the multi-component gaussian fit. For the errors in the peak 
intensities, T%, we consider the noise r.m.s (lcr rms ) of the 
spectra in Figure[T] and for the upper limits, the 3ov ms noise 
level of this emission. 



2.1 The broad SiO line emission toward 
condensations N, E and S 

Figure [Qshows that the brightest SiO J=2— >l and J=3— >2 
emission in the IRDC arises from three major condensa- 
tions toward (0,20) [N], (30,-30) [E] and (0,-80) [S]. The 
N and S condensations peak at the dense st cores reported 
in the filament (MM7/H 6 and MM6/H5; iRathborne et all 
12006] ; iButlerfe Tan|[2009l ), and harbor not only several 8 /im 
and 24 /xm sources, but 4. 5 fj,m extended emission (lik ely re- 
lated to H2 shocked gas; iNoriega-Crespo et all hoo4 ). This 
suggests that these condensations are active sites of star for- 
mation. 

Condensation E, however, is located east of the IRDC 
and slightly off the high extinction region (see Figure[T]) . Al- 
though small amounts of gas are pres ent, no local maxim a 
is seen toward E in either extinction dButler fc Tan! [2009). 
or in other high-density molecular tracers such as N2H + or 
H 13 CO+ (Caselli et al. 2010, in preparation). Like N and S, 
several 8 /im sources are detected toward this condensation, 
but only one is seen at 24 /im and this is about 15" north of 
the peak SiO position (see Figure[TJ. The Spitzer IRAC (3- 
8 fim) and MIPS 24 /jm fluxes of this source are consistent 
with a protostellar model with a luminosity of ~2x10 4 Lq, 
negligible circumstellar material, and a foreground extinc- 
tion of Ay=15 mag. This corresponds to a 15 Mq star on the 
zero age main sequence (ZAMS). This is the best fit model 



return ed from the SED fitting program of Robitaille et al.l 
(2007). However, this result is not unique and the deter- 
mination of protostellar properties from a relatively poorly 
constrained SED is quite uncertain. For example, the ob- 
served luminosity across the IRAC bands is only — 20 Lq, 
and the above estimate of a much higher luminosity relies 
on the reality of the large foreground extinction. The source 
in condensation E could be responsible for the SiO broad fea- 
ture observed toward this condensation. However, one would 
not expect strong outflow activity from such an evolved star. 
Alternatively, this source could be of much lower luminosity 
and mass and has more active accretion and outflow activ- 
ity. There may also be other low-luminosity protostars in the 
vicinity, although the lack of high enough angular resolution 
in our SiO images prevents to establish if this is the origin 
of the broad SiO emission in this condensation (see below) . 

In Figure[5] we show the P — V diagrams of the SiO 
J=2— >1 emission observed between the N and E, and the 
E and S condensations, for the velocity range from 40 to 
50kms _1 . Outside this velocity range, the emission of SiO 
.1=2— >1 is below the 2a rms level in the spectra. From 
Figure[5] we find that the typical SiO line profiles toward 
condensations N, E and S have a central component peaking 
at u_ts_R~45kms~ 1 , with broader line wing emission. The 
linewidths of the central component are ~l-3kms _1 , and 
those of the broad SiO emission are ~ 4-7 km s" 1 (Table©. 
The blue- and red-shifted terminal velocities are —40 and 
50kms _1 , respectively (i.e. ±5kms _1 with respect to the 
central velocity vlsr 1 ^ 45kms~ 1 ; Figure[T}. While the broad 
SiO emission toward condensation E is blue-shifted, the SiO 
linewings toward N and S mainly appear at red-shifted ve- 
locities. Broad line profiles have previously been reported in 
SiO and other molecular species such as H2CO and CH3OH 
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Figure 2. P - V diagrams of the SiO, 13 CO and C ls O J=2-»l 
emission (in units of T^), obtained along the lines linking con- 
densations N, E and S (see red dashed lines and red numbers in 
Figure[TJl. The distance (in arcseconds) shown in the j/-axis cor- 
responds to that measured along these lines from north to south 
(see also horizontal red dotted lines in this Figure). The contour 
levels for SiO J=2— >1 are 60 (2<r; dashed contours), 90, 120 and 
150 mK. For 13 CO and C 18 0, the first contour and level step are 
0.48 (3<t) and 0.8 K for 13 CO J=2-»l, and 0.48 (3a) and 0.48 K 
for C ls O J=2^1. While the beam size of the SiO J=2-s-l data 
is 28", the angular resolution used in the P — V diagrams of the 
13 CO and C ls O J=2->1 lines, is 11". 



toward several samples of IRDCs (see 



Carey et all 1 19981: 



20071 ; ISakai et al.1 



Beuther fc Sridharanl 120071 ; iLeurini et al. 
2008), and are believed to trace material associated with 
molecular outflows. 

The SiO P — V diagram of Figure[2] also shows that 
condensations N, E and S have several local maxima cen- 
tered at different radial velocities. Due to the low angular 
resolution of our SiO observations, it is currently impossi- 
ble to determine whether these maxima are produced by 
error fluctuations in the SiO line temperature, or whether 
they are associated with shocked gas high-velocity bullets, 
or with different low-mass protostars (see Section!?} . 

The P-V diagrams of 13 CO J=2-*T (Figure^ show 
that the terminal velocities of the 13 CO lines are, in gen- 
eral, similar to those measured for SiO toward the conden- 
sations N, E and S. The C 18 J=2— s>l emission does not 
reveal any significant broad line wing emission since it is 
mainly associated with the high-density gas seen in extinc- 
tion toward G035. 39-00. 33 (see Section [!0)) . The correlation 
betwee n the CO molecular g as and the mass surface density 
map of iButler fc Tan! (|2009l ). will be analysed in detail in 
the near future (Hernandez et al. 2010, in preparation). 



2.2 The extended and narrow SiO components 

In addition to the broad SiO condensation N, E and S, the 
P — V diagrams of Figure[2]show that very narrow SiO emis- 
sion [with linewidths of ^lkms' 1 ; offset (10,10) in Figure[T] 
arises from regions linking these condensations. This is more 
clearly seen in Figure|3] (central panel), where the narrow 
SiO lines arising from ambient gas (at 77Ls.R~45kms~ 1 ), 
form a large-scale and extended ridge that follows the fil- 
ament. We note that the SiO emission associated with 
the ridge shows narrower line profiles than C ls O J=2— !>1 
[0.8kms _1 vs. ~l-3kms _1 ; offset (10,10) in Table[2]. The 
peak intensity of narrow SiO toward this position is rela- 
tively weak (0.06 K; Table[5} and is at the ~5cr level (the 
a rma of the SiO J=2— »1 spectrum is 0.012 K; Figure[T}. 

The high-sensitivity single-pointing SiO spectra ob- 
tained with EMIR toward (-10,70) and (-37,37), also reveal 
that the narrow SiO component spreads north and northwest 
of Core MM7/H6 (Figure©. The SiO data toward (-10,70) 
has been smoothed to a velocity resolution of 0.53 kms -1 
to improve the signal-to-noise ratio of the spectrum. For 
the (-37,37) offset, however, we keep a velocity resolution of 
0.26 kms -1 , because the SiO line emission toward this po- 
sition is a factor of 2 narrower than that reported toward 
(-10,70) (i.e. 0.8 kms" 1 vs. 1.5 kms" 1 ; seeTable[2}. The nar- 
row SiO lines are very faint and have integrated line in- 
tensities of 0.038±0.005 Kkms" 1 and 0.030±0.005 Kkms" 1 
toward (-10,70) and (-37,37), respectively. S ince the SiO 
lines reported by iBeuther fc S ridharanl (|2007l ) in a sample 
of IRDCs have significantly larger linewidths (^2.5 kms" 1 ) 
than those observed toward the ridge, toward (-10,70) or 
toward (-37,37), these lines are the narrowest features de- 
tected so far in a high-mass star forming region. The SiO 
J— 3— >2 lines are not detected toward positions (-10,70) and 
(-37,37). 

From the weak intensity and spatial distribution of nar- 
row SiO, one could think that this emission could be due 
to some line emission contribution, within the large 30 m 
beam of our observations (~28"), arising from condensa- 
tions N, E and S. However, the narrow SiO lines measured 
toward (10,10), (-10,70) and (-37,37) have line profiles differ- 
ent from those observed toward N, E, and S (see Figure[T}, 
and their central radial velocities differ from those found to- 
ward these condensations. This is particularly clear toward 
(-10,70), where the line peak velocity of narrow SiO is red- 
shifted by ~1 kms" 1 with respect to that derived toward N. 
Therefore, the narrow SiO lines detected in G035. 39-00. 33 
trace different molecular material from that seen in the SiO 
condensations N, E and S. 

Finally, in Figure[T] we note that the narrow SiO emis- 
sion toward (-10,70) and (-37,37) lies below the 2a contour 
level of the integrated intensity SiO map. This is due to the 
fact that narrow SiO lines are diluted in the broad veloc- 
ity range considered to create the map. In Figure[3] (central 
panel), the narrow SiO emission detected toward (-10,70) 
with the high-sensitivity EMIR receivers (Figure[4]), also lies 
below the 3cr noise level of the SiO map at ambient veloc- 
ities, because the observations with the old SIS receivers 
were not sensitive enough to detect such faint emission (see 
above). The narrow SiO lines toward (-10,70) are of partic- 
ular interest because they do not show any clear association 
with a 8 jj,m or a 24 /im source. The narrow SiO emission 



6 I. Jimenez-Serra et al. 



detected toward (-37,37) is likely associated with the faint 
SiO condensation (intensity level of ~3<r) located at (-40,50) 
in the central panel of Figure[3] 

2.3 The CO filaments in G035. 39-00. 33 

The general kinematics of the low-density CO gas toward 
G035.39-00.33, as traced by C 18 J=2->1, are also shown 
in Figure[3] The CO gas is distributed along three differ- 
ent filaments with radial velocities «ls.r=44.1, 45.3 and 
46.6kms _1 . These values correspond to the averaged cen- 
tral radial velocities derived for every filament from the 
C ls O J—2~ >1 gaussian fit parameters shown in Tabled 
While the central and brightest filament at «lsa~45 km s _1 
bends east tightly following the densest material within 
the IRDC, the blue-shifted filament with wlsa~44 km s _1 
intersects the former one in an arc-like structure point- 
ing west (i.e. with the centre of curvature lying to the 
east of the IRDC). The two intersecting regions are co- 
incident with the highest density cores reported in the 
IRDC, MM7/H6, MM6/H 5 and MM8/H4 (|Rathborne et all 
120061 ; iButler fc Tad 12009 ). Although not as clear as for the 
blue-shifted CO component, the red-shifted filament with 
ii.ls.r~47 km s _1 seems to further bend east with respect 
to the central one, following the trend already shown by 
the latter filament. The morphology of these filaments re- 
sembles the molecular structures predicted by flow-driven, 
shock- induced and magnetically-reg ulated models of cloud 
form a tion (see e.g. iFiege fc Pudritz 2000a . bi ; Ivan Loo et al.l 
l2007i ; iHennebelle et al.ll200Sl ; iHeitsch etafll2009\ As dis- 
cussed in Sectional in these scenarios the high-density fila- 
ment seen in extinction would have been formed in a cloud- 
cloud collision, with the low-density CO filaments still ac- 
creting material onto it. The kinematics and excitation of 
the low-density CO gas will be described in a forthcoming 
paper (Jimenez-Serra et al. 2010, in preparation). 



3 EXCITATION AND FRACTIONAL 
ABUNDANCES OF SIO 

Except for the (10,10) position, where SiO J=3— >2 single- 
dish observations were not carried out (see Section[2]), we 
can derive the H2 densities of the gas, n(H2), and the ex- 
citation temperatures, T ex , and column densities of SiO, 
iV(SiO), from the SiO J=3->2/ J=2->1 intensity ratio and 
by using the Large Velocity Gradient (LVG) approximation. 
As shown in Section HOl the validity of this approximation 
is provided by the fact that the SiO line emission toward 
G035. 39-00. 33 is optically thin. We assume that the SiO 
J=2— >1 and J=3— >2 lines trace the same material. 

3.1 LVG modelling and input parameters 

In our LVG calculations, we have considered the first 15 ro- 
tational levels of SiO, an d used the col li sional coefficients of 
SiO with H2 derived by iTurner et al.l (| 19921 ) for tempera- 
tures up to T=300K. These collisional coefficients are well 
suited for our case of study, because the kinetic temperatures 
assumed for the SiO shocked gas in the IRDC are T^C300 K. 
As input parameters, we have considered the CMB (Cosmic 
Microwave Background) temperature (i.e. T(, 9 =2.7K), the 



linewidths of the SiO emission, and the kinetic temperature 
of the gas, Tkin- The brightness temperature, Tb, excita- 
tion temperature, T ex , and optical depth, r(SiO), of every 
rotational line transition of SiO, are then calculated for a 
grid of models with different H2 (volumetric) gas densities 
and column densities of SiO. We note that T ex may signif- 
icantly differ from Tkin, since the excitation conditions of 
the SiO shocked gas could be far away from the LTE (Local 
Thermodynamic Equilibrium). 

Three different velocity regimes have been considered 
for the SiO emission observed in G035. 39-00. 33: the ambi- 
ent component, which ranges from 44.5 to 45.5kms _1 ; the 
moderate shocked gas, with 42.5 kms _1 ^ULg_R<44.5 kms -1 
or 45.5 kms~ 1 <«LSi{^47.5 kms -1 ; and the high- velocity 
regime, with «ls_r<42.5 kms -1 or i>ls.r>47.5 kms -1 . 
Within every velocity regime, we perform the LVG cal- 
culations for velocity bins of lkms -1 . We assume that 
Tfcm = 15K for the ambient comp onent (simi l ar to those 
derived in a sample of IRDCs; iPillai et al l I2007T ). and 
Tfcm=25K and Tfe 4n =45K for the moderate and high- 
velocity SiO gas, respectively. The latter two temperatures 
are consistent with those found in the shock ed gas of the 
L1448-mm outflow (| Jimenez-Serra et al]|2005h . 

We would like to stress that the selection of these tem- 
peratures are expected not to be crucial in our calculations 
of JV(SiO), since the excitation of the SiO rotational lines 
with J up <5 does n ot strongly depend on Tkin but on the H2 
density of the gas jNisini et al.l |2 007 i). Indeed, if we increase 
Tkin from 15 K to 50 K for the narrow SiO component, and 
from 25-45 K to 300 K for the moderate and high- velocity 
SiO shocked gas, the derived iV(SiO) change by less than a 
factor of 2. 

The brightness temperatures, Tb, derived with the LVG 
model for the SiO J=2— 5*1 and J— 3— >2 lines, are finally com- 
pared with those observed toward G035. 39-00. 33, in units 
of T m b- We assume that the SiO emission is uniformly dis- 
tributed and that the beam-filling factor of our SiO obser- 
vations is ~1 (only in this case Ts^Tnb)- This is justified 
by the fact that the SiO emission is extended in the IRDC. 

To derive the SiO fractional abundances, the H2 col- 
umn densities were estimated from C ls O for the ambient 
gas, from 13 CO for the moderate velocity regime, and from 
CO for the high- velo city gas. In con t rast w ith the mass sur- 
face density map of IButler fc Tanl (|2009h , which gives an 
averaged value of the H2 column density of the gas toward 
a certain position, CO and its isotopologues provide esti- 
mates of the H2 column density as a function of velocity 
within the shock. We assume isotopic ra tios 12 C/ 13 C=53 
and 16 0/ ls O=327 (|Wilson fc Roodlll994h . and a CO frac- 
tional abundance of 2xl0~ 4 across the IRDC. The uncer- 
tainty in the CO abundance is about a factor of two, consid- 
ering its variations in different molecu lar cloud complexes 
(e.g. lFrerkingT"L angcr. & Wilson 1982 ) as well as w ithin 
the same complex (jPineda. Caselli fc Goodmanll2008r i. The 
value adopted in this study is close to that directly measured 
towards an other high mass star forming region (NGC2024; 
lLacv et al]|l994l 'l. which better represents the properties of 
IRDCs. 
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Figure 3. Integrated intensity maps (in units of kms -1 ) of the SiO J=2— >1 (black contours) and C ls O J=2— >1 lines (red contours) 
measured toward G035. 39-00. 33 from 40 to 44.5kms _1 (blue-shifted gas; left panel ), 44,5 to 45.5 kms - 1 (ambient gas; central panel), 
and 45.5 to 50 kms -1 (red-shifted gas; right panel). The mass surface density map (Butler & Tan 2009) is shown in greyscale (contour 
levels as in FigurefTJ- The contour levels of the SiO J=2— >1 maps are 24 (3a), 32, 40, 48 and 56 mKkms for the blue-shifted emission, 
51 (3<r), 68, 85, 102, 119 and 136 mKkms" 1 for the ambient velocity range, and 27 (3a), 36, 45, 54, 63, 72 and 81 mKkms- 1 for the 
red-shifted gas. For clarity, we only plot the contours at the 5a, 15a and 25<r levels in the C ls O maps, with first contours at 150, 350 
and 200mK kms -1 for the blue-shifted (left), ambient (centre), and red-shifted (right) velocity ranges, respectively. Beam sizes of the 
SiO J=2— >■! (black circle) and C ls O J=2— »1 observations (filled red circle) are indicated in the upper left corner. 



3.2 LVG results 

The T ex and SiO fractional abundances, x(SiO), derived to- 
ward (-10,70), (-37,37), (0,20), (30,-30) and (0,-80) by means 
of the LVG approximation, are shown in Figure^ Toward 
(10,10), x(SiO) was calculated assuming that T ex =9K [i.e. 
similar to those derived toward (-10,70) and (-37,37); see 
below] and that the SiO emission is optically thin. 

For the narrow SiO emission toward ( — 10,70) and 
(-37,37), we obtain T ex ~9K, optical depths r(SiO)<0.01, 
and SiO column densities ranging from 5xl0 10 cm- 2 
to 10 11 cm -2 . This implies SiO fractional abundances 
x(SiO)~6-7xlO -11 (FigureEJ). The derived H 2 gas densi- 
ties are ^6xl0 5 cm -3 . 
fractional abundance is 



Toward (10,10), the derived SiO 



-5x10" 



These abundances are 



a factor of 10 la rger than the upper limits found in dark 
cloud s (^10~ 12 ; IZiurvs et al.l 11989 ; iRecmena- Torres et al.l 
2007), and are similar to those measured from narrow SiO 
toward the molecular outflows in the low-mass star form- 



ing regions NGC1333 an d L1448-mm (|Lefloch et al.lll998l ; 
Ijimenez^ Serraet al.|[20ol ). 

Toward the N, E and S condensations, T ex and x(SiO) 
tend to increase from the ambient to the moderate and the 
high- velocity regimes (Figure^). The typical optical depths 
derived for the SiO emission in these condensations are 
r(SiO)^0.06. The derived Eh gas densities and SiO col- 
umn densities range from 10 5 cm -3 to 10 6 cm -3 , and from 
5xl0 10 cm~ 3 to 4xl0 11 cm~ 3 , respectively. For T ex , a sim- 
ilar behavior for the excitation of the SiO shocked gas has 
been reported toward th e L1157-mm and L1448-mm out- 
flows jNisini et al.ll2007t ). where the SiO J=8^7/J=5^4 
ratio is known to increase as a function of velocity within the 
shock (|Nisini et al.l 12007). In the case of the SiO fractional 
abundances, x(SiO) is progressively enhanced from ~10~ 10 
in the ambient gas, to ~10~ 9 in the moderate velocity com- 
ponent, and to JSlO -8 in the high- velocity shocked gas. This 
trend has also been observ ed toward the L1448-mm outflow 
ijjimenez-Serra et al"1l2005l ). 
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Although the LVG model does not provide the errors 
associated with T ex , we can roughly estimate them from 
the integrated intensity ratio between the SiO J=3— >2 and 
the SiO J=2— >1 transitions. The errors range from 15% to 
35%. This range corresponds to the change in signal-to- noise 
ratio going from the strong SiO lines at ambient velocities to 
the fainter SiO emission at the moderate and high-velocity 
shock regimes. 
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Figure 4. High-sensitivity SiO J=Z—¥\ spectra measured with 
EMIR toward (-10,70) and (-37,37). The rms is 6mK for the 
(—10,70) spectrum and 10 mK for the (—37,37) observation. 



3.3 Validity of the LVG approximation 



As shown in Section [3T2l the typical optical depths derived 
for the SiO emission toward G035. 39-00. 33, are r(SiO)«l. 
In molecular outflows, the SiO line emission is expected to be 
optically thin for ratios A'(SiO)/Au<5x 10 13 cm" 2 km" 1 s, 
where TV(SiO) is the derived SiO c olumn density an d At), 
the linewidth of the SiO line profiles (|Nisini et al.ll2007n . The 
derived SiO column densities, A r (SiO), toward the IRDC are 
relatively low and range from 5xl0 10 cm" 2 to 4xlO n cm" 2 
(Section l3.2|) . Considering that the velocity bins used in 
our calculations are lkms _1 -wide, the ratio iV(SiO)/Au is 
5xl0 10 -4xl0 1 1 cm" 2 km" 1 s, i.e. well below the upper limit 
established bv lNisini et al.l l|2007l ). Therefore, the use of the 
LVG approximation in our case is fully justified. 



3.4 



Comparison of the SiO fractional abundances 
with shock modelling 



By us ing the parametrized model of Ijimenez-Serra et al.l 
(2008) that mimics the steady-state physical structure of a 
perpendicular C-shock, we find that the low SiO fractional 
abundances of ^lO" 11 associated with the narrow SiO emis- 
sion in G035. 39-00. 33, could be generated by the sputtering 
of dust grains within C-shocks with ^ s ~12kms _1 . A mini- 
mum ion-neutral drift velocity of v c i=\v n — i>j|=5kms is 
required to sputter SiO abundances of a few 10" 11 from the 
icy mantles of dust grains, if a fraction al abundance of 10~ 8 
for S iO is present within the mantles l|jimenez-Serra et al.l 
2008). This SiO abundance is the same as that measured in 
the moderate velocity SiO shocked gas in L1448-mm, and 
believed to trace m aterial recently released fr om the man- 
tles of dust grains (| Jimenez-Serra et al.l 120051 ') . If SiO were 
only present within the cores of dust grains, a shock ve- 
locity of Us^Okms" 1 , with an ion- neutral drift velocity of 
v,i~18kms" 1 , would be required to reproduce the low SiO 
fractional abundances of few 10" 11 for narrow SiO. 

The progressive enhancement of the SiO fractional 
abundances from ~10~ 10 to Jsl0~ 8 toward the N, E and 
S condensations, could be explained by the sputtering of 
the mantles of dust grains along the propagation of a C- 
shock w ith i> a ~30kms~ , as prop osed for the L1448-mm 
outflow (I Jimenez-Serra et al.l 120081 ). However, we note that 
other mechanisms suc h as shattering or vaporis ation in non- 
dissociative J-shocks (|Guillet et al.ll2007l , l2009l ), could addi- 
tionally contribute to the processing of dust grains in these 
regions. The inclusion of those mechanisms in the currently 
available state-of-the-art shock descriptions, might signifi- 
cantly change our understanding of the shocked SiO emis- 
sion in molecular outflows. 

The comparison of synthetic SiO line profiles (as derived 
from radiative transfer modelling) with the SiO line emission 



observed toward the IRDC, lies outside the scope of the 
present paper. The lack of high enough angular resolution 
in our SiO OTF maps, indeed prevents to clearly establish 
the actual morphology and size of the SiO emitting regions, 
which are crucial parameters in the modelling of the SiO 
line profiles. 



4 ON THE ORIGIN OF THE PARSEC-SCALE 
SIO EMISSION IN G035.39-00.33 

Theo ret ical models o f flow-d riv en llHennebelle fc Perauit] 



19991; iHeitsch et all 120061 ; iHennebelle et all 



2008 



Heitsch et al.l 120091). shock-induced ( Kovama fc Inutsukal 



200d. 120021 ; Ivan Loo et "ail 120071), and magnetica lf 



-irt" 

regulated formation of clouds IjFiege fc Pudritzl l2000al lbl) , 
predict that these regions have a very filamentary structure 
at their early stages of evolution. Consistent with this idea, 
the filamentary IRDC G035. 39-00. 33 shows a relatively high 
number of quiescent cores (without H2 shocked gas or 24 (im 
sources), w hich are believed to b e at a pre-stellar/cluster 
core phase (|Chambers et al.l 12009). As a consequence, one 
should not expect to find a significant impact of outflow 
inte raction within the cores and on their surroundin gs (see 
e.g. iMartm-Pintado et al.lll992l ; iBeuther et~alll2002l ). 

The high-sensitivity maps of the SiO emission toward 
G035. 39-00. 33, however, reveal for the first time that SiO 
is widespread along an IRDC. Large-scale SiO emission 
(with sizes ranging from 4 to 20 pc) has also been re- 
orted across the molecular clouds in the Galactic Ce ntre 
Martm-Pintado et al.lll997l ; lAmo-Baladron et al.ll2009l ). In 
this case, the origin of these lines is different from that 
in G035. 39-00. 33, because the SiO gas is highly turbu- 
lent (linewidths of ~60-90 kms" 1 ). The large SiO fractional 
abundances (~10~ 9 ) derived toward these regions are likely 
generated in fast shocks of superno va explosions, HII re- 
gions and Wolf-Rayet stellar winds |Martm-Pintado et al.l 
Il997l ). and /or associated with X-ray or cosmic ray induced 
chemistry (|Amo-Baladr6n et al.|[20CM ). Besides the GC, the 
large-scale SiO emission observed toward this IRDC consti- 
tutes the largest SiO feature detected so far in a star forming 
region. 

In Section[2j we have shown that the SiO line profiles 
measured toward G035. 39-00. 33 have two different compo- 
nents with different spatial distributions, kinematics and ex- 
citation. The first one consists of bright and compact SiO 
condensations (N, E and S) with broad line profiles in both 
SiO J=2— >1 and J=3— »2 transitions. From our excitation 
and fractional abundance analysis of the SiO lines, T ex 
and x(SiO) tend to progressively increase for larger veloci- 
ties within the shock (from 6 to 40 K, and from ~10~ 10 to 
>10" 8 , respectively), as expected for shocked gas in molec- 
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Figure 5. Excitation temperatures, T ex (color lines), and SiO 
fractional abundances, x(SiO) (symbols), derived for the ambi- 
ent, moderate and high-velocity regimes (vertical dotted lines) 
observed toward several offsets in G035. 39-00. 33. We consider 
«0=45.0kms — 1 . The errors associated with T ex are estimated 
to range from 15% to 35%. 



ular outflows (|jimenez-Serra et al ] |2005l : iNisini et alJ [2007) . 
Toward N and S, this idea is s upported by the detection 
of 8//m sources, 24 ^m sources JCarev et al. I 2001SI ) and ex- 
tended H2 shocked gas ((Chambers et al.ll2009h closely asso- 
ciated with these condensations. In addition, the observed 
spatial extent of the broad SiO emission is only ^30" (i.e. 
^0.4 pc at 2.9 kpc), which is consistent with those mea- 
sured in outflows associated with high-m ass protostars to- 
ward other IRDCs (|Fallscheer et aljfeoogl ). We also note that 
this size is similar to the beam size of the SiO J=2— sT 
observations, for which the fractional abundances of SiO 
have been reproduced by our C-shock and sputtering models 
(Sectionl3^l. 

For condensation E, we have seen that the broad 
SiO emission is not associated with any massive core in 
the region, and is clearly off the high e xtinction region 
ijRathborne et alJl200d : IButler fc Tad" l2009). In addition, no 
local maxima is seen in extincti on or in other high -density 
molecular tracers in this region jButler fc Tani 2009, Caselli 
et al. 2010, in preparation). However, the IRAC 3-8 /im and 
the MIPS 24 /im fluxes of the source located 15" north of the 
SiO peak are consistent with a 15 Mq star in the ZAMS. Al- 
though one would not expect strong outflow activity from 
such an evolved star, this source could be the origin of the 
broad SiO emission found in this condensation (Section l2.ip . 
Indeed, the SiO fractional abundances derived toward E are 
the largest measured in G035. 39-00. 33, and are consistent 
with the idea that the SiO gas has been processed by a rel- 
atively strong shock (with ii s ~30kms _1 ; Section l3~4]> . sim- 
ilar to those found in outflows associated with high-mas s 
star formation (see e.g. IRAS 20126+4104; ISu et al.ll2007l ). 
Therefore, although with the current data it is not possible 
to clearly establish the origin of broad SiO toward conden- 
sation E, it is likely that this emission is related to shocked 
outflowing gas. 

The second SiO component detected toward G035.39- 
00.33, consists of extended and very narrow SiO emis- 
sion that links condensations N, E and S, and which also 
spreads very faintly toward the north and northwest of 
Core MM7/H6. From the low SiO fractional abundances de- 
rived for this component (~10 -11 ), the narrow SiO emission 
could be produced by the interaction of a large-scale low- 
velocity shock with ii s ~12kms _1 (Section I3.4[) generated 



boost iHeitsch et al.ll2009T ). In G035.39-00.33. this collision 
could have been produced by the interaction between the 
main filament (as s een in the mass surface density map of 
IButler fc T an 2009), and the lower density filaments traced 
by C ls O (see Figure[3]). Since the time-scales required for 
SiO to freeze-out onto dust grains are relatively short (i.e. 
from 5 x 10 3 yr to 5 x 10 4 yr for volume densities from 1 5 to 
10 4 cm" 3 ; see Section 6.2 in lMartfn-Pintado et aHll992T ), the 
dust grain processing event associated with this interaction 
would be relatively recent. The narrow feature would then 
constitute a signature of the filament-filament collision or 
of previous accretion events produced onto the main IRDC 
filament. 

This scenario is supported not only by the extended 
morphology of narrow SiO, but by the fact that the co- 
herent CO filaments observed toward G035. 39-00. 33 re- 
semble the mo l ecula r structures predicted in these models 
van Loo et al. 120071 ; iHennebelle al.l 2008; He itsch et al.l 



l2009h . In addition, the narrow SiO lines, specially toward 
the north and northeast of Core MM7/H6, do not show any 
clear association with 8 /j,m or 24 /im sources fSection l2.2[l . 



However, we note that the relatively small (2-4kms _1 ) ve- 
locity difference between the CO filaments compared to the 
shock velocity of ~12kms _1 required to produce the low 
SiO fractional abundances of ~10 -11 , does place constraints 
on this scenario. This would require either that much of the 
relative velocity between the colliding molecular gas is in the 
plane of the sky or that much of the gas has already been de- 
celerated in the interaction. Detailed comparison with the 
results of numerical simulations are required to assess the 
likelihood of these possibilities. 

As proposed bv lLefloch et~ail (|l998T > for the NGC1333 
low-mass star forming region, the narrow SiO emission in 
G035. 39-00. 33 could arise from decelerated shocked gas as- 
sociated with large-scale outflows driven by the 8 /jm and 
24 /im sources seen in the IRDC. This gas would have been 
decelerated by its i nteraction with a de nse and clumpy sur- 
rounding medium |Lefloch et al.| [l998). It is also possible 
that narrow SiO is produced by material recently processed 
in the magnetic precurso r of MHD shocks, as pro posed 
for the L1448-mm outflow ( Jimenc z-Serra et al. 2004) . This 
idea i s similar to that suggested by Beuther fc Sridharanl 
(2007), for which narrow SiO would be linked to the 
youngest jet/outflow objects present in their sample of 
IRDCs. As discussed by these authors, narrow SiO lines 
are unlikely to be produced by an effect of outflow incli- 
nation with respect to the line-of-sight, because this would 
lead to the detection of fewer outflows with broad line 
emission than outflows wit h narrow SiO (see Section 3.2 in 
iBeuther fc Sridharanl 12007] ) . Although any of the previous 
mechanisms could explain the narrow SiO lines in the ridge 
between condensations N, E and S, it seems unlikely for off- 
sets (-10,70) and (-37,37), where the narrow SiO lines do 
not show a clear association with 8 /im or 24 /jm sources 
(Section[23|). 

Alternatively, the extended and narrow SiO emission 
toward G035. 39-00. 33 could be produced by a widespread 
and lower-mass population of protostars, compared to those 
powering condensations N, E and S. Some of these dis- 
tributed protostars may be visible as the 8 fim sources in 



10 /. Jimenez-Serra et al. 



Figure[T] although we do see SiO emission from regions ap- 
parently devoid of such sources. In this scenario, beam di- 
lution would then prevent us from detecting the broad SiO 
line wings expected to arise from these objects. Interfero- 
metric observations are thus needed to discriminate between 
i) the large-scale shock scenario, remnant of the IRDC for- 
mation process; ii) decelerated or recently processed gas in 
the precursor of MHD shocks in large-scale outflows driven 
by the 8 fim and 24 /xm sources; and iii) an undetected and 
widespread lower mass protostar population, as an origin of 
the widespread SiO emission in G035. 39-00. 33. 

Extended narrow SiO emission could also be produced 
by the UV photo-evaporation of the mantles of dust grains 
in photon dominate d regions (PDRs) such as the Orion bar 
l|Schilke et al.ll200lf) . However, this mechanism seems un- 
likely in G035. 39-00. 33, because of the relatively low lumi- 
nosity of the region (it is observed as an infrared dark rather 
than bright cloud), and because the UV radiation field re- 
quired to produce SiO fractional abundances similar to those 
observed in this cloud (of ~10 -11 ), sh ould be at least few 
hundred times the Galactic UV field l|Schilke et al.ll200ll) . 
There is no evidence for sources capable of producing such 
an intense FUV field in this region, and even propagation of 
the much lower-intensity Galactic FUV radiation field into 
this cloud, would be strongly impeded by its high extinc- 
tion. Nor are cosmic ray induced UV photons l|Gredel et al.l 
Il989h expected to play a key role in the formation of narrow 
SiO since the same UV field is generated in nearby quiescent 
dark clouds, where no SiO is detected. 

In summary, we report the detection of widespread (size 
of ^2pc) SiO emission toward a very filamentary IRDC. 
This emission presents two different components with dif- 
ferent kinematics, excitation and spatial distributions. The 
compact morphology, large SiO fractional abundances, and 
broad SiO line profiles observed toward N, E and S, indicate 
that these condensations are shocked gas in outflows associ- 
ated with high-mass star formation. The second SiO compo- 
nent is extended along the filament and shows very narrow 
line profiles, low SiO abundances, and lower excitation than 
the gas detected toward N, E and S. Although interfero- 
metric images are needed to clearly establish the origin of 
this emission, the properties of narrow SiO are consistent 
with i) a large-scale shock, remnant of the IRDC formation 
processes; ii) decelerated or recently shocked material in the 
precursor of shocks in large-scale outflows powered by 8 /im 
and 24 /im sources; or iii) an undetected and widespread 
population of lower mass protostars. 
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